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ABSTRACT 

We report unsuccessful searches for pulsations from the neutron star RX J1836. 2+5925 identified 
with the EGRET source 3EG J1835+5918. A 24-hr observation with the NRAO Green Bank Telescope 
at 820 MHz placed an upper limit on flux density of 17 /iJy for P > 10 ms, and gradually increasing 
limits for 1 ^ ^ 10 ms. The equivalent luminosity is lower than that of any known pulsar with the 
possible exception of the radio-quiet 7-ray pulsar Geminga. A set of observations with the Chandra 
X-ray Observatory HRC totaling 118 ks revealed no pulsar with 1 ms < P < 10 s. The upper limit 
on its pulsed fraction is 35% assuming a sinusoidal pulse shape. The position of RX J1836. 2-1-5925 in 
Chandra observations separated by 3 years is unchanged within errors, leading to an upper limit on 
its proper motion of < 0.14" yr~^, or vt < 530 km s~^ at d = 800 pc, a maximum distance estimated 
from its thermal X-ray spectrum. With these null results, the properties of 3EG J1835-I-5918 and its 
X-ray counterpart RX J1836. 2-1-5925 are consistent with a more distant or older version of Geminga, 
or perhaps a recycled pulsar. Having nearly exhausted the capabilities of current instrumentation at 
all wavelengths, it will likely fall to the Gamma-ray Large Area Space Telescope to discover pulsations 
from 3EG J1835+5918. 

Subject headings: gamma rays: observations — stars: neutron — X-rays: individual 
(RX J1836.2-f 5925) 



1. INTRODUCTION 

The brightest of the " unidentified" high-G alactic- 
latitude EGRET sources (Hartman^efaD 119991), 3EG 
J1835-f5918 at {£,b) = (89°, +25°) has long been 
associated with the X-ray emitting neutron-star 
RX J1836. 2+5925, which remains its most plausible 
counterpart even though pulsations have not been de- 
tected at any wavelength. Unlike blazars, which are 
highly variable and have steeper 7-ray spectra, 3EG 
J1835+5918 shows no evidence for long-term variabil- 
ity, and its spectrum can be fitted by a relatively flat 
power law of photon index F = 1.7 fro m 70 MeV to 
4 Ge V, with a turn-down above 4 GeV (Reime r et al.l 
120011) ■ similar to known 7-ray pulsars. RX J1836. 2+5925 
is the only optically undetected X-ray source in the 
7-ray error box, all t he others being cla ssifled as un- 
likely 7-ray emitters (|Mirabal et al.|[2000f) . The detec- 
tion of RX J1836. 2+5925 as a weak, ultrasoft source 
in th e ROSAT All-Sky Survey (,Mirabal fc Halp ern 
|2001| ) suggested that it is a thermally emitting neu- 
tron star that is either older or more distant than 
the 7-ray pulsar Geming a ([Halpern fc RudermanI 119931 : 
iBignarni fc Caraveol[T99l . 

Observations using the Chandra Advanced CCD 
Imaging Camera (ACIS) and the Hubble Space Tele- 
scope (HST) furth er supported this interpretation 
([Halpern et al.l[2002f l. Two components were required 
to fit the Chandra X-ray spectrum, a blackbody of T 
3.5 X 10^ K and a power law of photon index F — 2.2+0.6. 
The non-thermal X-ray component is important evidence 
of magnetospheric activity that is found in all 7-ray pul- 
sars, but not in all cooling neutron stars. An optical 
upper limit of V > 28.5 from the HST verified the 
neutron-star interpretation of RX J1836. 2+5925 from its 
extreme X-ray-to-optical fiux ratio, fx/fv > 6000. We 



also used the thermal X-ray spectrum and optical limit 
to bound the distance to RX J1836. 2+5925 in the range 
250 < d < 800 pc.y The energetics of 3EG J1835+5918 
are plausible for a pulsar at d < 800 pc since its 7-ray 
luminosity (assumed isotropic) is 3.8 X 1034(d/800pc)2 
ergs s~^, comparable to the spin-down power E = —Iflfl 
of Geminga (3.3 x 10^"* ergs s~^). Efficiencies approach- 
ing 100% are achieved by those 7-ray pulsars having the 
s mallest spin- down pow er . 

iHalpern et al.l ([20021 ) also searched for radio pulsa- 
tions from RX J1836.2+5925 using the JodreU Bank 
Lovell telescope at a frequency of 1.4 GHz, achiev- 
ing an upper limit of S'1.4 < 0.1 mJy. Adopting 
an upper limit of d < 800 pc, the pulsed pseudo- 

< 



luminosity limit of RX J1836. 2+5925 is L1.4 = SiAd 
0.064 mJy kpc^. There a re only four pulsa rs known 
with L1.4 < 0.1 r nJy kpc^ (ICamilo et al.ll2002[ ) in addi- 
tion to Geminga ([McLaughlin et anil999L and references 
therein). 

In this paper, we report an even more sensitive search 
for radio pulsations from RX J1836. 2+5925 using the 
NRAO Green Bank Telescope (GBT), as well as the first 
X-ray pulsation search with the Chandra High Resolu- 
tion Camera (HRC-S). These investigations approach the 
sensitivity limits of existing instrumentation to search for 
the pulsar in 3EG J1835+5918. 

2. GBT RADIO PULSAR SEARCH 

RX J1836. 2+5925 is circumpolar at the GBT, where 
on 2002 December 6 we observed it for 24 hr using the 
BCPM spectrometer (i Backer et "allll997[ ). Due to the in- 
termediate latitude of the target, with a relatively small 
predicted dispersion measure (DM) and low Galactic 
background temperature, we observed at a central fre- 
quency of 820 MHz, with a bandwidth of 48 MHz in each 
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TABLE 1 

Chandra HRC Observations of RX J1836. 2+5925 



ObsID Start Time Live Time Counts 
OJT^ 

6182 2005 Feb 28 19:27 45,018 284 

4606 2005 Mar 9 15:30 28,030 189 

6183 2005 Mar 11 08:19 44,999 256 

Total 118,047 729 



of two polarizations. The signals from corresponding po- 
larization channels were summed in hardware, and the 
total-power samples from each of 96 frequency channels 
were sampled every 72 fis and recorded to disk for off-line 
analysis. 

We analyzed the data using standard pulsar search 
tech niques, implemented in the PRESTO software pack- 
age ([Ransom et al.l [200^ . First we identified and ex- 
cised radio-frequency interference. We then dedispersed 
the data in the DM range 0-110 cm~^ pc, with resolu- 
tion of 0.3 cm^'' pc. This is twice the m aximum Galactic 
DM p redicted for this line of sight by the lCordes fc Laziol 
(|2002f ) electron density model; for the distance range of 
250-800 pc, within which the neutron star likely lies (§1), 
the predicted DM range is 2-9 cm~^ pc. Each time se- 
ries contained 1.2 billion points, and would have been 
needlessly challenging to analyze. We therefore down- 
sampled the data to a resolution of 0.288 ms, resulting 
in much more reasonable time series of 300 million points. 
For a DM of 40 cm^'^ pc, the smearing caused by disper- 
sion across one of the frequency channels is 0.3 ms, so 
that for smaller DMs the effective time resolution of the 
search was about 0.3 ms. 

We performed the periodicity search using fast-Fourier 
techniques, and were sensitive to periods P ^ 1 ms. In 
order to improve sensitivity for a range of pulse duty cy- 
cles, we incoherently summed up to 16 harmonics of the 
data; for P > 10 ms, we maintained good sensitivity for 
duty cycles down to about 0.03P. The sensitivity at long 
periods decreases gradually due to red noise, but this is 
not a significant concern for P < 5 s. Unfortunately, we 
found no convincing pulsar candidates in this search. 

Using the standard modification to the radiometer 
equation, assuming a pulse shape with 10% duty cy- 
cle, we were sensitive to pulsars with P > 10 ms hav- 
ing a period-averaged flux density at 820 MHz > 17 ^Jy, 
with gradually decreasing sensitivity for smaller peri- 
ods. Converting this to the more usual 1.4 GHz pulsar 
search frequency wi th a typical spectral index of -1.6 
(jLorimer et al.|[l995l ) gives Sia < 7 /iJy. For a distance 
of 800 pc, this corresponds L1.4 < 5 iiJy kpc^ . All kno wn 
radio pulsars have a greater L1.4 (jCamilo et al.l[200^ . 

Based on these results, it is therefore unlikely that 
RX J1836. 2-1-5925 is a radio pulsar beaming toward the 
Earth. More likely, it is either beaming away, or is not 
an active radio pulsar at all. 

3. CHANDRA HRC PULSAR SEARCH 

Three observations of RX J1836.2-H5925 totaling 
118 ks of live time were performed between 2005 Febru- 
ary 28 and 2005 March 11 using the Chandra HRC-S in 
timing mode, which has time resolution of 16 /iS, and 
provides absolute accuracy of 42 /is as determined from 
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Fig. 1. — The parameter space {f,f) searched for X-ray pul- 
sations from RX J1836. 2+5925 is shaded. The step size 5f = 
4 X 10"^"^ Hz s~^, and the search includes / = 0. Lines of con- 
stant E and Tc are shown. The star indicates the parameters of 
Geminga. 

observations of the Crab pulsar^ . Table [1] lists the ba- 
sic parameters of the observations, including the number 
of photons extracted from a l."2 radius aperture. The 
live time is 99.6% of the elapsed time, limited by event 
processing. The total number of photons detected from 
RX J1836. 2-1-5925 is 729, of which 20 are estimated to 
be background (using nearby, source-free regions of the 
image). The count rates from the three observations are 
within 9% or less of the average count rate, 0.006 s~^, 
consistent with no variation. The spatial distribution of 
photons is consistent with a point source. Originally con- 
ceived as a contiguous, 135 ks pointing, the observation 
could not be executed as planned because of new oper- 
ational restrictions that came into effect. The resulting 
large gaps in time diminish the sensitivity to pulsations 
because the increased number of independent frequencies 
and frequency derivatives that have to be searched yield 
more false candidates. 

The total time span of the observations from beginning 
to end is T = 9.6 x 10^ s, which determines the inde- 
pendent intervals for frequency and frequency derivative, 
A/ = 1/2T and A/ = l/T^, respectively. We chose a 
frequency step size Sf = 2 x 10~^ Hz to oversample the 
independent frequency interval by a factor of « 2.5, and 
(5/ = 4 X 10""'^^ Hz to achieve a similar oversampling 
factor for frequency derivative. 

We transformed the photon arrival times to Barycen- 
tric Dynamical Time, and did coherent searches of 
the full data set using the discrete Fourier transform 
technique, also kno wn as the Rayleigh or Zf test 
(jBuccheri et al.|[r983f ). We searched the entire range of 
(/, /) parameter space for a pulsar that has 0.1 < / < 
1000 Hz, characteristic age Tc = —f/2f > 20 kyr, as 
well as £^ = —An'^Iff < 10'^'' ergs s^^. In practice, 
we found it simplest to confine the search to the region 
shaded in Figure [U bounded by/< 1.4x 10~^^ Hzs~^ 
for / < 17.9 Hz, and / < 2.5 x 10-^° f'^ Hz s^^ for 
/ > 17.9 Hz. The total number of trials thus defined is 

http://cxc.harvard.edu/cal/ 
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« 1.6 X 10 , an oversampling of the independent trials 
by a factor of « 7. 

This is a liberal search range that is justified as follows: 
RX J1836. 2+5925 would have a blackbody temperature 
greater than its observed ^ 3.5 x 10^ K if it were younger 
than 20 kyr; an independent argument for an older age 
is given in §4. It would have a wind nebula or stronger 
nonthermal X-ray component if > 10'^^ ergs s~^. In 
view of the distance limit d < 800 pc, the ratio L^/E < 
10"^ if _E = 10'^^ ergs a~^, a smaller ratio than all other 
pulsars, which typically have lO""' < L^/E < 10~^. So 
it is unlikely that E is as large as 10^^ ergs s~^. 

The range of parameters searched also includes all of 
the known millisecond (recycled) pulsars, and is effective 
as long as the pulsar does not have a neutron star binary 
companion. The deep HST limit, corresponding to ab- 
solute magnitude My > 19 at d = 800 pc, rules out the 
more common white dwarf companions. 

No significant pulsed signal was detected. The largest 
values of Zf found were « 44. The theoretical dis- 
tribution of follows that of with 2n degrees of 
freedom. For n = 1, the distribution is an exponen- 
tial with a mean of 2, so the single-trial probability that 
Zf > 44 by chance is 2.8 x 10~^°. This is expected to 
arise r andomly in « 10^° independent trials. iLeahv et al.l 
([1983') showed that to detect sinusoidal pulsation with 
amplitude (pulsed fraction) fp at a power level S = Z^ 
with 50% probability, the number of photons needed 
is iV = 25'//2. In our case, N = 729 (709 corrected 
for background), so we find fp < 0.35 corresponding to 
Zi < 44. We verified this analytic expression by exam- 
ining folded light curves at the periods associated with 
the maximum values of Z^. 

It is not surprising that the pulsed fraction of a 7- 
ray pulsar in X-rays should be less than 35%. This 
is not a very sensitive limit consdering that other 7- 
ray pulsars whose soft X-rays are predominantly ther- 
mal (Geminga, PSR B1055-52, possibly PSR B0656+14) 
have even s maller pulsed fractions (30%, 21%, and 12%, 
respectivelv. iDe Luca et all l2005l ) . In these cases, the 
pulse profiles are quasi-sinusoidal. The Z\ test is there- 
fore a good one for RX J1836.2-f 5925, whose X-rays are 
also dominated by a soft thermal source within the re- 
sponse of the Chandra HRC, although the result is not 
very restrictive in this case. We also note that the Vela 
pulsar has a more compli cated pulse shape, while its 
pulsed fraction is only 9% ((Manzah et al.ll2007f l. 

4. CHANDRA ASTROMETRY 

In order to test for proper motion of the neutron 
star, we compared the position of RX J1836. 2-1-5925 
on the HRC with the ACIS position obtained 3 years 
earlier (jHalpern et al.l I2002D . First, we updated the 
ACIS position to the USNO-Bl.O system using opti- 
cally identified sources in the ACIS image. The result is 
R.A. = 18''36'"13."685, decl. = -H59°25'29."95 (J2000.0), 
and differs from the previously used USNO-A2.0 sys- 
tem by 0."3. Then we registered the HRC image to the 
ACIS image using seven X-ray sources in the vicinity of 
RX J1836.2-f5925, which required a zero-point shift of 
the HRC image of +0.''052 in R.A. and -t-0."89 in deck. 
The resulting HRC position of RX J1836. 2-^5925 is R.A. 
= 18'^36'"13.''674, decl. = -|-59°25'30."15 (J2000.0). 



The difference between the positions of 
RX J1836.2+5925 in 2002 and 2005 is therefore 
only 0."21, which is comparable to their combined 
statistical errors. In order to bound the possible proper 
motion, we adopt an upper limit that is twice this 
difference, or 0."42, corresponding to /i < 0."14 yr~^. 
This allows the neutron star to have traveled from a 
birth in the Galactic plane to its present position in 
6.4 X 10^ yr or longer, which is not an unreasonable age 
for its X-ray temperature. At the maximum distance 
of 800 pc, the limit on /i corresponds to a tangential 
velocity limit vt < 530 km s~^, which is typical for radio 
pulsars. 

5. DISCUSSION AND CONCLUSIONS 

Even in the absence of detectable pulsations, 
RX J1836. 2-1-5925 remains the leading (and only) candi- 
date for identification with 3EG J1835-h5918. Its soft X- 
ray spectrum, and absence of optical and radio emission, 
support the hypothesis that it is an older and possibly 
more distant cousin of the Geminga pulsar. The upper 
limit on proper motion that we derived is a significant 
new observational constraint, allowing an age of at least 
6.4 X 10^ yr for RX J1836.2-h5925 ff it was born in the 
Galactic plane. This compares favorably with the char- 
acteristic age Tc = 3.4 x 10^ yr of Geminga, and is con- 
sistent with a surface temperature of only ^ 3.5 x 10^ K, 
compared with 4.8 x 10^ K for Geminga (ICara veo et al.l 
120041 : iJackson fc Halperi3l2005t iKargaltsev et a l. 2005). 

Assuming that the age and distance of 
RX J1836.2+5925 are both larger than those of 
Geminga, it is likely that 3EG J1835-I-5918 is a max- 
imally efficient 7-ray pu lsar close to its death line 
(jChen fc Rudermaiilll993D . following the trend of 7-ray 
efficiency increa sing with decreasing spin-down power 
(iThompson et a l. 1997, 1999). Alternatively, we cannot 
rule out that RX J1836. 2-1-5925 is a millisecond pulsar 
with similar spin-down power and magnetospheric gap 
voltage as Geminga. If so, its thermal X-ray emission 
may be due to surface reheating by a magnetospheric 
accelerator. 

6. PROSPECTS FOR GLAST 

It is not now possible to make a significantly more sen- 
sitive search for X-ray pulsations from RX J1836. 2+5925, 
as the operational limitations of Chandra no longer allow 
long, contiguous pointings. While XMM-Newton is in 
principle more sensitive, it has not been used to observe 
RX J1836. 2+5925, which lies near the Earth-avoidance 
zone of the satellite orbit, with maximum allowed point- 
ing durations of « 20 ks. 

Finally, the upcoming Camma-ray Large Area Space 
Telescope (CLAST) may be the best hope of detecting 
pulsations from 3EG J1835+5918. GLAST wiW certainly 
reduce the positional uncertainty of 3EG J1835+5918, 
which will be an independent test of association with 
RX J1836. 2+5925, and additional motivation for a dedi- 
cated pulsar search. For this, a s ufficiently long, pointed 
observation is best, as argued bv iRanso a ([200^ . 

We estimate the expected CLAST count rate from 
3EG J1835+5918 using its photon fiux of 5.9 x 10"'^ 
cm~^ s"^ between 100 MeV and 4 GeV (jReimer et al.l 
|2001|) . Assuming an average on-axis effective area of 
5000 cm^ for the CLAST Large Area Telescope (LAT) in 
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this energy range, the predicted count rate is 0.003 s . 
A more detailed calculation, convol ving the spect r um o f 
3EG J1835+5918 from Figure 3 of iReimer etlHI (|200lD 
with the LAT on- axis effective area curve^, predicts a 
similar rate of 0.0027 s'^ in the 100-4000 MeV range. 
This is almost half the count rate from the Chandra 
HRC, and holds great promise for a pulsar discovery. 

It is likely that GLAST will succeed where Chandra 
failed because, unlike thermal X-rays, 7-rays have sharp 
peaks with pulsed amplitude approaching 100%. In 1 
week of elapsed time, including Earth blockage, CLAST 
can collect as many photons from 3EG J1835-I-5918 as 
we have obtained using Chandra. On-source time in a 
pointed observation of 3EG J1835+5918 can be higher 
than average because the source is close to the ecliptic 
pole. Also, Galactic background is negligible at these co- 
ordinates {i,b) = (89°, +25°). If GLAST is successful. 



the Chandra photons analyzed here may yet be used to 
find a coincident period and period derivative with confi- 
dence, thus establishing the identity of 3EG J1835-I-5918 
with RX J1836. 2+5925, and determining its physical pa- 
rameters. 
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